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Abstract 

The evolution equations of Einstein's theory and of Maxwell's theory — the latter 
used as a simple model to illustrate the former — are written in gauge covariant first 
order symmetric hyperbolic form with only physically natural characteristic direc- 
tions and speeds for the dynamical variables. Quantities representing gauge degrees 
of freedom [the spatial shift vector x^) and the spatial scalar potential (j){t, x^), 
respectively] are not among the dynamical variables: the gauge and the physical 
quantities in the evolution equations are effectively decoupled. For example, the 
gauge quantities could be obtained as functions of {t, x^) from subsidiary equations 
that are not part of the evolution equations. Propagation of certain ("radiative") 
dynamical variables along the physical light cone is gauge invariant while the re- 
maining dynamical variables are dragged along the axes orthogonal to the spacelike 
time slices by the propagating variables. We obtain these results by (1) taking a 
further time derivative of the equation of motion of the canonical momentum, and 
(2) adding a covariant spatial derivative of the momentum constraints of general 
relativity (Lagrange multiplier /?*) or of the Gauss's law constraint of electromag- 
netism (Lagrange multiplier (j)). General relativity also requires a harmonic time 
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slicing condition or a specific generalization of it that brings in the Hamiltonian 
constraint when we pass to first order symmetric form. The dynamically propa^ 
gating gravity fields straightforwardly determine the "electric" or "tidal" parts of 
the Riemann tensor. 

PACS #: 04.20.-q, 97.60.Lf 
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1 Introduction 



We examine the Cauchy problem for general relativity as the time history of the 
two fundamental forms of the geometry of a spacelike hypersurface, its metric g and its 
extrinsic curvature K. By using a 3 + 1 decomposition of the Riemann and Ricci tensors, 
we split the Einstein equations into initial data constraints — equations containing only 
g, K, and their space derivatives — and evolution equations giving the time derivatives 
of g and K in terms of space derivatives of these quantities and also of the lapse and 
shift, that is, of the variables that fix the proper time separation between leaves of the 
foliation of spacetime and the "time lines" threading the foliation 0, ^, Q . The constraints 
can be posed and solved as an elliptic system by known methods. (See, for example 
[|I], |, 1^.) However, the equations of evolution of g and K are in essence the spatially 
covariant canonical Hamiltonian equations with first order time derivatives and second 
order space derivatives. Casting these equations into hyperbolic form |]^, ^, |^, ^ is 
a problem important for theoretical analysis and for practical applications. However, 
the procedure used may entail a loss of spatial coordinate covariance, the use of non- 
geometrical variables, or the appearance of unphysical characteristics, that is, directions 
which are neither along the light cone nor orthogonal to the time slices. 

All of these difficulties can be overcome, as we show, by obtaining first a nonlinear 
wave equation for the evolution of K. The evolution of the spatial metric g is just its 
dragging by K along the axis orthogonal to the spacelike time slices. 

One method for obtaining this result requires choosing the time slicing by specifying 
the mean (extrinsic) curvature of the slices. It leads to a mixed elliptic-hyperbolic system 
for which we have recently proven local in time, global in space, existence theorems in 



the cases of compact [O] or asymptotically fiat slices [13, 12, 14 



In this paper, we concentrate on a second method, which relies on a harmonic time 



slicing condition or a certain class of generalizations of it |T^, |T^, |Tj, |T6[. We obtain 
equations of motion equivalent to a covariant first order symmetric hyperbolic system 
with only physical characteristics. We construct this system explicitly. The only quanti- 
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ties propagated along the light cone are curvatures. The space coordinates and the shift 
vector are arbitrary. In this sense, the system is gauge covariant. 



2 Electromagnetism 

Maxwell's theory provides a simple example for illustrating the method we use in treating 
Einstein's theory. The idea is to work with the dynamical or 3+1 form of the theory, that 
is, essentially the canonical form. By using a further time and a further space deriva- 
tive in a definite way, one proceeds to construct a physically natural second order wave 
equation for the canonical momentum, in this case the electric field. The wave equation 
for the electric field can be written in a first order symmetric hyperbolic form ( "FOSH" 
form) that is also flux conservative. Furthermore, this procedure maintains gauge covari- 
ance and produces a FOSH system that has only the physically natiu^al characteristic 
directions given by null generators of the light cone and timelike vectors orthogonal to 
the spacelike time slices. We refer to these properties as "simple physical character- 
istics." Correspondingly, one obtains "characteristic flelds" that propagate along the 
light cone. By treating this simple case flrst, many steps in dealing with the technically 
more intricate case of general relativity can be abbreviated because they follow the same 
pattern. 

Consider Maxwell's theory on flat spacehke shces of flat Minkowski spacetime with 



and prescribed current jn{t,x^) = {—p,ji) such that V^j^ = dop + V*ji = (an overbar 
denotes a spatial tensor or operator). Using F^,^ = 2V[^A,y] and and = {—(pjAi), we 
deflne i?^ = V^F^i^, which leads to the identities 



ds 
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dt^ + Qijdx' 



(1) 



(2) 



Ri = -doEi + V'Fij. 



(3) 
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The standard form of the Maxwell equations can be written as 

Ro + 47rp = 0, (4C) 

Ri - ATTji = 0, {4E) 

doAi = -Ei - (4L>) 

where Vj^ = c?j0. ( C denotes "constraint," E "evolution," and D "definition.") To 
construct the new 3+1 system, form the spatial one-form 

= doRi - ViRo (5) 

and use (4D) and the definition of F^j to obtain 

ni = {-dodo + A)Ei = aEi, (6) 

where A = V'^Vfc. While (6) is an identity, the Maxwell equations give us a new state- 
ment of the evolution equation, 

Q, = 47r(9oj- + V,p). (7) 

The new form of the Maxwell system is therefore 

V'Ei = 47rp, (8C) 

nEi^A7r{doji + Wip), (SW) 

doAi = -Ei - Vi(t>, {AD) 

where W denotes "wave equation," and the scalar potential = x^) is an arbitrary 
gauge field that can be regarded for the present as prescribed. 

The initial data (Cauchy data) for the new system on a spacelike slice Sq (t = 0) 
are (a) given freely, (b) E^ such that (8C) holds for given p, and (c) d^Ei such that 
(4E) holds {Ri — Anji — 0) for given jV In other terms, the extra time derivative used in 
constructing requires that we give as further Cauchy data doEi satisfying on Eq the 
usual Maxwell equation of motion, which is (3) combined with (4E). 



For any given (j){t,x^), the new system (8W), (4D) forms a quasi-diagonal strictly 
hyperbolic system in the unknowns Ei and Ai. Its characteristics at a point are the light 
cone and the time axis do, which lies inside the cone. The cone therefore determines 
the propagation. The new system in its present form implies that propagates at light 
speed along the cone while is dragged along the time axis. 

Before passing to FOSH form, let us show that the new system is equivalent to the 
standard Maxwell equations: Assume (7) has been solved. Then 

V% = 47rnp, (8) 

where current conservation, V^j^ = 0, has been used. Therefore, using (6) we find 

□ (V^Ei -47rp) = 0. (9) 

But V'-Bj — 47rp = {Rq + 47rp = 0) on Eq by hypothesis, and from the "Bianchi 
identity" V^-R^ = (i.e. V^VF^,, = 0) and current conservation Vj^ = 0, we have 

doiV'Ei - 47rp) = V'{-Ri + 47rj-) = (10) 

on Eq. In the last equality we have used the initial data condition (c) that Ri — Airji = 
on Eq. Therefore, the Cauchy data pertaining to (10) are zero and it follows that the 
Maxwell constraint (4C) holds on Eq x I, where I is an interval of time. (Mathematically 
rigorous statements of results in this paper can be given in terms of Sobolev spaces Hg, 
or weighted Sobolev spaces Hg^s for asymptotically flat time slices.) 

Now (5) and (7) together with the result just obtained that the constraint (4C) holds 
on Eq X I show that 

9o(i?,-47rj-) = (11) 

on Eo X I. But again R^ — iirji = on Eq, so it must remain zero on Eq x I, which is 
what we wanted to prove. 

In passing to the FOSH form of the Maxwell equations, we will leave the constraint 
(8C) unchanged; it can be converted to an eUiptic equation that needs to be solved only 
on Eq. Also, the other initial-value data will be given as above. 
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We consider here the vacuum case p = ji = 0. Defining Goi = d^Ei and Gji = VjEi, 
we write the first order system corresponding to (4D) and (8W) as 

doAi = -E, - {AD) 

doEi = Goi, {13D) 

^oGo^ - V'G,, = 0, il3W) 

doGji - VjGoi = 0, (13A) 

where A= "adjoint", an equation obtained by using the trivial commutator [So, Vj] = 0. 
(It is aheady clear that this system has "simple physical characteristics" as previously 
defined.) 

Define the transpose of a column matrix u by 

= (A, E, Go, Gi, G2, G3) , (14) 

where the final index "i" has been omitted for each entry. The above equations (4D)- 
(13 A) have the form 

dou + B^VjU + C = 0, (15) 

where B-' (j = 1, 2, 3) is an ri x n square matrix and C is an n x 1 column vector, with 
n the number of unknowns. (With the abbreviated notation (14), n = 6.) The matrices 
B-' and vector C depend on {t,x^) in the linear case and on (u;t,x*) in the quasi-linear 
case. Equations (15) are hyperbolic and symmetrizable if there is an n x n symmetric 
positive definite matrix M'' such that M'^B-' = M-' is symmetric. (It may be that these 
properties hold only in some restricted region of spacetime and, in the quasi-linear case, 
only in a functional neighborhood of a particular u). If we carry out this operation, we 
obtain an explicitly FOSH form, 

M"V<,u + M°C = M°9ou + M-' VjU + M°C = 0. (16) 

From the FOSH form, one can construct "energy norms" and "energy inequalities" that 
lead, by the methods of functional analysis, to precise results on the existence and unique- 



ness of solutions, and global estimates of the time of existence |]TT 



In the present case, the B-^'s are not symmetric but the symmetrizer is 

„ / 1,3x3 03x3 \ , ^ 

M° = _ 1 (17) 

\O3x3 §3x3/ 

where (g'~^) = g^^ ■ The matrix M° is positive definite if and only if g^j is properly 
Riemannian. (This is important in the case of gravity, where gij is an unknown. There, 
the symmetrizer will be essentially two copies of centered on the diagonal of a simple 
larger matrix.) We verify easily that the M^B-^'s are symmetric. 

It is clear from the form of (4D)-(13A) that and are fields associated with 
evolution along the the time axis. There are also fields propagating at light speed. These 
"characteristic fields" can be found by returning to the symmetrizable form (15). Then, 
to exhibit one example of such fields explicitly, we choose a space direction at an event 
corresponding, for example, to j=l. This means the direction ±ni, orthogonal in space 
to the — constant two-surface passing through the event and lying in the t — constant 
slice containing the event. We solve the eigenvalue problem associated with B^ and 
diagonalize it with an appropriate matrix T by a similarity transformation 

TB^T^^ = D\ (18) 

where is diagonal. The non-zero eigenvalues of B^ give the speeds of the "charac- 
teristic fields" with respect to observers at rest in the time slices, whose four-velocities 
are parallel to the orthogonal time axis ("Eulerian observers"). By construction, for 
the choice of direction ±ni, these characterisitic fields include only those which propa- 
gate in the local spacetime two- plane defined by and ±ni at the given event. As we 
change the direction ±ni, we move around the locally isotropic null cone at the event 
and capture all the characteristic fields. 

The non-zero eigenvalues of B^ are ±q;~^, where a = (g^^)'^^"^ is the spatial "lapse 
function" corresponding to the chosen direction. The quantity a~^, of course, is the 
speed of light because, in this example, the time is proper but the spatial displacements 
are not. (In local Cartesian coordinates we obtain, of course, ±1.) 
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The characteristic fields are obtained from 

TU = Uehar. (19) 

They are found to be 



—;=a 
V2 



'G^:,±a-\Gu-v''Gai)\^ (20) 

where v"" = —g^"'{g^^)^^ with a = 2, 3 is the "shift vector" in space corresponding to the 
X =constant two-surface whose "space lapse function" is a = (g^^y^/^. Substituting 
into (20) the definitions of Goi and Gji, we see that the characteristic fields in this example 
are the differential electric field along the advanced and retarded "±ni- branches" of the 
generators of the null cones. That is, the fields (20) are 

[do ± a~\V, - v'^Va))] E,. (21) 

The operator in rectangular brackets (times l/-\/2) gives the geometrical form of the 
usual advanced and retarded time derivatives defined with respect to the chosen local 
spacetime two-plane at the given event. Knowing the characteristic fields facilitates 



imposing physically consistent boundary conditions [T^, in this case, boundary conditions 
compatible with causality. As examples, consider a "no outgoing waves" condition at the 
future event horizon of a black hole or a "no incoming waves" condition in the vacuum 
far- field of electromagnetic source switched on a finite time in the past. 

The final points we wish to emphasize concern the facts that (1) no particular gauge 
condition was imposed on Ai to obtain our physical FOSH system, and (2) the scalar 
potential ("gauge field") Aq = —0 played an entirely passive role. It could either be given 
or obtained (say) by a supplementary (non-Maxwell) equation, as a function of (t, x*). Any 
mathematically consistent method of obtaining 0(t, x*) can be used without affecting the 
features of the Maxwell equations per se in the form in which we have given them above. 
The "gauge sector" and the "physical sector" have thus been separated to this extent. 
Our FOSH system and its simple physical characteristics are maintained in the face of 
complete gauge covariance. (A similar statement holds for Ao in the Yang-Mills equations 
[IT5| and for the "shift- vector" of the Einstein equations — see below.) 



3 Gravitation 



In Einstein's theory it is convenient to write the metric as 

ds'' = -N\ey + g.jO'e^, (22) 
with 9^ = dt and 9^ = dx^ + (3'^dt, where is the shift vector. The cobasis 9°' satisfies 



d9'' = -^C^^9f^ A 9^ (23) 



with Cqj = — Cjg = djP^ and all other structure coefficients zero. The corresponding 
vector basis is given by Cq = dt — P^dj and Cj = di where dt = d/dt, di = d/dx\ and 
the action of Cq on space scalars is the Pfaffian or convective derivative, eo[/] = d^f = 
dtf - f3^d,f. 

We shall assume throughout that the lapse function > and the space metric g 
on St is properly Riemannian. Note that cjij = Qij and g^^ = g'^K 
The spacetime connection one-forms are defined by Q 

^ ■'"^7 + 9"^ '^sip9-f)e - ^C'^^ = + ^ip^], (24) 

where (P^y) = ^(/37 + 7/5), [P'j] = |(/37 — 7/3), and T denotes a Christoffel symbol. These 
connection forms are written out in [Q. In particular, the extrinsic curvature (second 
fundamental tensor) of the space slices is given by 

K,^ = -Nul = -^N-'^og^J, (25) 

where we define for any t— dependent space tensor T another such tensor d^T of the 
same type by setting 

4 = ^ - Cp, (26) 

where is the Lie derivative on St with respect to the vector /3. Note that do and di 
commute: [do, di]T = 0. 

We have chosen in this paper to work with Kij rather than the usual canonical 
momentum, (det gY^'^{Hg^^ — K^^), conjugate to gij, where H = K\ denotes the mean 
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curvature. This is a matter of choice, not of principle. Indeed, our general procedure 
can be applied to other choices of canonical or geometric variables. (We note that Kij is 
conjugate to {detgY^'^g'^K) 

We regard the dynamics of general relativity as being given by the time dependence 
of the first and second fundamental forms of the foliation S^. It is therefore appropriate 
to use (9o as our time derivative and Vi, the spatial covariant derivative, for spatial 
derivatives = r*;, = f*;;.). An advantage of using do is that it is always timelike 
because it is orthogonal to the spacelike slices S^. We recall that d/dt is not necessarily 
timelike in Einstein's theory so that using it to define hyperbolicity can be confusing. 
Unlike the flat-spacetime Maxwell theory on flat, spacelike hypersurfaces, here we have 
the non-trivial commutator 

[4, V,]m, = Nuk [2V(,ir,)^ - V^K,, + 2a(,ir,)' - a^K,^ , (27) 

where Uj is a spatial one-form and aj = di log = Vj log is the acceleration one-form 
of the "Eulerian" observers at rest in the time slices. The commutator on higher-rank 
tensors is similar to (27), with corresponding terms for each further index. 

It is necessary to relate the spacetime and space Riemann tensor components as well 
as the corresponding Ricci tensor components. We use the conventions in and the 
0" cobasis. The Riemann tensor components are given by 

Rijki = Rijki + '^Ki[kKiy, (28) 
Roijk = 2NVijKk]i, (29) 
RoiOj = N0oKi,+NKikK'',+W^WjN), (30) 



where we note that 



NViVjN = N^iViaj + aiaj), (31) 



and, in three dimensions, 

Rijki = '29i[kRi]j + '2gj[iRk]i + R9i[i9k]j- (32) 
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The latter relation proves to be very important in our construction of a FOSH system 
for the Einstein equations in (3+1) dimensions. 

The Ricci curvature Rap = aap is given in (3+1) form by the identities 

Rij = Rij- N-%Kij + HKij-2KikK''j-Viaj-aiaj, (33) 
i?oi = NV^{Hg,j-Kij), (34) 
Roo = NAN -N'^KmkK'^^ ^Nd^H. (35) 

The Einstein equations, Gap — i^Tap, can be written as Rap — upap, where k — Sir ii 
G — c — 1, and, in (3+1) dimensions, pap — Tap — \gapT'^a- In this paper we will 
take Tap = but indicate at certain stages what the matter terms are. (The presence 
of matter fields indeed causes no difficulties in a complete hyperbolic system for gravity 
plus matter if these fields themselves have well-posed Cauchy problems in fiat spacetime 
and, as a sufficient condition, are minimally coupled to gravity.) 

We recognize (33) as the evolution equation E, (34) as the "momentum constraints" 
Cj, and from the definition Gap = Rap — \Rgap of the Einstein tensor, we see that the 
"Hamiltonian constraint" Cq is defined by 

G\ = i(i?% - R\) = \{K^kK^'' -H^- R). (36) 

The constraint equations Cj and Co can be posed and solved as an elliptic problem and 
will not be discussed further here. Therefore, in analogy to (4C), (4E), and (4D), we 
have 

i?oi = 0, ^iR%-R\)^0, (37C) 

Rij = 0, (37^;) 
doQij = -2NKij. (37D) 

To construct the new (3+1) system, form the symmetric spatial tensor 

^ij = doRij — 2V(ji?j)o, (38) 
12 



and work it out using the relations above and 

^oR^, = Vk{do^)-VjidoTi), (39) 
4fJ,. = -2\/^.iNK,)') + V\NK,,), (40) 

to obtain, in analogy to (6), the identity 

Qij = NuKi, + J^j + Si,, (41) 



where □ = —{N ^d^Y + This wave operator was denoted "□" in [|1^; but here we 
reserve that standard symbol for (^^^VaV/j, where is the usual spacetime covariant 



derivative operator. In |]14| and ||T2[, slightly different wave operators were used, resulting 
in slightly different counterparts for Jij. These wave operators all have the same principal 
symbol. 

With the conventions above, to which we adhere in this paper, the explicit form of 



the " source" Jij is given in ||15|| ; below, it will be written out explicitly for our FOSH 
system. For the present, we note that Jij contains no derivatives of Kij higher than the 
first; and that it contains Rijki and Rij with the former subject to elimination in favor 
of the latter using (32). The contribution of matter to (41), before the elimination of 
Rij to be carried out later in obtaining FOSH form, is found by noting that Einstein's 
equations assert 

VLij = K[doPij - 2V(iPj)o]. (42) 
The "slicing term" Sij is given by 

Si, = -N-'V,V,{doN + N^H). (43) 

Because of the presence of the VjVjiJ term in Sij, the principal part of the complete 
second-order operator acting on Kij in (41) is not a wave operator. However, this problem 
(which has no counterpart in Maxwell's theory) can be eliminated by using either the 

doN + N^H = 0, (44) 
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simple harmonic slicing condition {g°'^uj^ij = 0) 



or a generalized harmonic slicing condition |]T5 



doN + N^H = Nf{t, x^), (45) 



where f{t, x^) is an arbitrary given function analogous to a "gauge source" Although 
we shall, for simplicity, take / = in most of the following discussion, a choice of / 7^ 
causes no material changes. In any event, we now have, from (41) and (45), 

= NuKij + Jij - ViVjf - 2a(iV,)/ - /(V^a, + a,a,), (46) 

which is a second-order nonlinear wave operator on Kij for all time slicings in the class 
(45). (Recall at = VilogN.) 

[Another way to handle Sij is to prescribe the mean curvature: set H = h{t,x), a 
known function. This leads to a mixed elliptic-hyperbolic representation of the Einstein 
evolution equations which we will not pursue here. See 0, |12|.] 



The new Einstein system of equations of motion is given by (46), (45), and dogij = 
—2NKij. The appropriate Cauchy data on an initial slice Sq (t = 0) are g, K, N, and 
SqK. We require g and K to satisfy the usual momentum and Hamiltonian constraints 
and OqK. to be given by Rij = in (33) in the vacuum case: that is, we require that all 
the Einstein equations G^fs = hold on Sq. The lapse function must satisfy > 0. Note 
that the shift which we may regard at present as being known in the form /5*(t, x-'), 
is playing a passive role completely analogous to the scalar potential (p in our treatment 
of the Maxwell equations. The shift has been absorbed into the operator do- 

To show that the new system is indeed solvable, before showing that it is equivalent to 
the Einstein equations on Sq x /, we shall reduce it to the form (8W) and (4D) obtained 
for the Maxwell equations. To accomplish this, we use the expression for H in terms of 
gij and its derivatives [cf. (37D)] to write (45) in the form 

do\og[N{detgr'/'] = f{t,x^), (47) 

from which A^ can be obtained in terms of the spatial metric and the known functions 
and /, called "algebraic gauge" in |1^, (Compare [0] in the case = and / = 
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and [0, [T^ in the case / = 0.) The new system now acquires a reduced form 



oKi, = -N-'[J,, ^ViVjf -2a^S^j)f - fiy^aj + aia,)], (48) 
do9ij = -2NKij, (49) 

where N is understood to be determined by (47). These equations, analogous to the 
Maxwell equations (8W) and (4D), are a quasi-diagonal system whose characteristics are 
the light cone and the do axis, which lies inside the cone. This form of the equations 
implies that Kij propagates along the light cone at unit speed while gij is dragged by Kij 
along the Oq axis. Since equivalence to Einstein's theory has not yet been proven and 



therefore cannot be used, the Ricci curvatures occuring in Jij ||T5[ are understood by their 
definition in terms of the metric gij and its derivatives, up to second order. Using (49) in 
(48), we obtain a third order strictly hyperbolic equation for {det g)~^^'^ gij. A theorem 



of Leray|T9| then implies the temporally local existence of solutions in spatially local 



Sobolev spaces. In practice, we see that we can regard (45), (46), and dogij = —2NKij 
system to be solved simultaneously. 

It is next required that we show the equivalence of the above system to the standard 
Einstein equations. This proof is transparent especially when the actions of the operators 
do and Vi on spatial scalars, one-forms, and covariant symmetric tensors are related to 
those of Vo and Vi on these same objects. Such relations follow from the definitions of 
the operators and the values of the connection one-forms u^^. 

Assume that the reduced system defined by (45) or (47), as well as (48) and (49), has 
been solved, so that the lapse and the metric are known on So x I, where I is an interval 
of time. We must prove that given the Cauchy data for the reduced system, i.e., that the 
Einstein equations hold initially and that > is given initially in accord with (47), 
then (38) [and (42)], supplemented by the twice-contracted Bianchi identities [and the 
covariant conservation law V^Tajj = 0], implies that they hold on Sq x I. It is necessary 
that we have already found (i.e. proven the existence) of the lapse and the metric so 
that (38) becomes a definite equation for the Ricci curvature, and we can raise and lower 
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indices as necessary. In vacuum, 

^ij = doRij - 2V(,i?,)o = WoRij - 2W^iGj)o + (l.h.E.)o = 0, (50) 

where the notation (l.h.E.)„ denotes additive terms that are hnear and homogeneous in 
the Einstein tensor and its derivatives of order < n. Form the tensor 

= Qij - QijQ^ = VoG,,- - 2V(,G',)o + gijV'Gko + (l.h.E.)o. (51) 
From the twice-contracted Bianchi identities V^G^i, = 0, we have 

VOGoo + V^G.o = 0, (52) 
V^Goi + V'Gji = 0. (53) 

Compute W^fl[j, use the Ricci identity to commute covariant derivatives, and use (53) to 
obtain 

□G'oi + (l.h.E.)i-0. (54) 

But Gai3 = on Eo by hypothesis and dgGa/s — (or VoGq^ = 0) on Eq from (51), 
(52) and (53). Hence Goi = on Eq x I. (The action of do on Gap is defined because in 
our basis, Gqo is a space scalar, Goi is a space one- form, and Gij a space tensor.) 

From Goi = on Eo x I and (51), it now follows that doGij = 0; therefore, Gij = on 
Eo X 1 because it was zero on Eq. Likewise, from (52), we infer that Goo = on Eq x I, 
which completes the proof. 

The passage to FOSH form of the Einstein equations proceeds similarly to that for the 
Maxwell equations except for two points. One is technical: the non-trivial commutator 
(27) has to be used. The other is conceptual: the necessity of using a generalized harmonic 
time-slicing condition. In what follows, we shall use simple harmonic slicing (44) in the 
vacuum case; use of (45) only adds known terms and no new variables. The key point 
about harmonic slicing doN + N'^H = is that when it is combined with the Einstein 
equation it!oo = 0, i.e. 

doH = -AN + K^kK'^'^N, (55) 
16 



we obtain a wave equation for A^: 

□A^ - {KijK'^ - H^)N = 0. (56) 

It is at this stage that the Hamiltonian constraint enters exphcitly into our system: 
see (36). By applying Vfc to this equation we obtain a non-hnear wave equation for a^. 
We may then reduce the pair of wave equations involving nKij and na^ each to flux- 
conservative FOSH form as we did previously for OEi. Furthermore, the new variables we 
must introduce to carry this out are precisely the ones required to express the non-linear 
right-hand-sides of these two wave equations in terms of the defined variables, given the 
elimination of Rij in these equations using the standard Einstein equation of motion (33) 
with Rij = 0. {Rijki is to be regarded as eliminated from the right-hand-side of the 
nKij equation using (32): this key step in obtaining FOSH form uses the fact of three 
space dimensions. The nKij equation itself has the same form in any Lorentzian (d+1) 



spacetime for the Einstein equations.) These steps are executed in [|1^, |T5|, |T^, [T^, but 
here we shall state the final results slightly more explicitly. 

The new variables needed to reduce the nKij and naj wave equations are, respectively, 

Loij = N'^doKij Lkij = VkKij, (57) 

aoi = N'^d^ai aji = VjQi. (58) 

(In previous works we denoted L^ij by Lij and L^ij by M^ij.) 

Altogether, our unknowns up to this point in the analysis are N, gij, Kij, ai, aoi, aji, 
Loij, and Lkij. Note that the shift is not among the unknown variables. Like (j){t,x^) 
in the case of Maxwell's equations, f3\t,x^) may be regarded as given or obtained (as a 
function of t and x^) in any manner that is mathematically compatible with the Einstein 
equations in FOSH form. Obtaining (3^{t,x^) is needed in solving Einstein's equations, 
but its "manufacture" is not a part of Einstein's equations per se. 

The reduction to first order form produces the following system (in vacuum, with 
simple harmonic slicing) 

d^N = -N'^H, (59) 
17 



Oij 



doQij = -2NKij, (60) 

doKij = NLoij, (61) 

doLoij - NV^Lkij = NJij, (62) 

+2K^(, {K^)"'ak + a^)K^^ - a^K,)^)] , (63) 

dotti = Naoi, (64) 

+a, + a^ak + 2a\ - 2KmkK^'')] , (65) 



dottij - NVjaoi = N 



(66) 



ttk {2L(^ij)^ — L^ij + 2a{iKj)^ — a'^Kij^ + a^a^i 

where, in (62), NJ^j — J^j with insertion of the unknowns and ehmination of spatial 
Riemann and Ricci tensors: 



+ak [SL'^ij - 4L(ij)'=) + 4a(iL^)/ - (lOLofe(i + WKkmK'^a + lOaka^i + 8ak(^ Kj)^ (67) 
+2gijK"''' (^Lornk - HK^k + 2KmiK\ + a^ak + amk) - QijH {K^kK""^ - H"^^ . 

Note that the right-hand sides of (59)- (66) are polynomials in the unknowns and in 
g^K However, the reduction to first order form is not quite complete because the spatial 
covariant derivatives on the left-hand sides of (62), (63), (65), and (66) involve derivatives 
of the unknown g^j. Two ways to complete the reduction are as follows. 

One could introduce a known fiducial torsion-free connection -Fjfe(^) ^0 on Eq xl whose 
associated covariant derivative operator is (say) Di. Then define the tensor 



sik^rik-Fikit.A 

as a new unknown whose equation of motion is 



(68) 



(69) 
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where 

doT]^ = -N[2L(,u) - L'',k + 2a^^Kk)' - a^Kj^]- (70) 

The equation obtained from combining (69) and (70) can be appended to (59)- (66) to 
form the complete first order system. The operators Vj are replaced by relations of the 
form 

VkTij = DkTij - TijSl^ - TiiSj^. (71) 

Another method that can be used to complete the reduction is simply to append the 
relation (70) and to rewrite the Vj's in terms of ordinary derivatives: all the equations 
of the system are still spatially covariant though some will have, of course, individual 
terms which are not tensors. This is true, perhaps not obviously, of (70) also. One has 
to recall that 

CfsT], = {p%^,-^,^lP^+fl,^JP'+^^^,p')+^J^,p' (72) 

is a spatial tensor. 

In our system, (62) and (65) are wave equations and (63) and (66) are, respectively, 
their "adjoint" equations. All the other equations are definitions except for the slicing 
equation (59) and the equation for the evolution of the connection. This system is 
symmetrizable upon multiplication by a symmetric positive definite matrix consisting of 
I's down the main diagonal and two sub-matrices of the form (17) centered on the main 
diagonal to handle the two wave equations and their adjoints. The unknowns can be 
expressed by 

= (g, K, Lo, Li, La, L3, N, a, ao, ai, aa, as, T), (73) 

where a trailing pair of indices is suppressed from g, K, and L^, a trailing vector index 
from a and a^ and all of the indices from T. Therefore, we have constructed a complete 
FOSH system. The structure of the left-hand sides has the "flux-conservative" form 
possible for curved space and curvilinear coordinates. 
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For completeness, we state explicitly the initial data for the FOSH system (59)- (66) 
and (70), given /5*(t, x^). These data are consistent with those of (47), (48) and (49) with 
/ = 0. The lapse function A'^ > is chosen freely, which implies the initial acceleration 
Qi — di log N. The metric gij and extrinsic curvature Kij satisfy the usual Gauss-Codazzi 
constraints C — and Ci — given by (37C). From gij and we obtain the initial 
values of f*-;;. and of Lijk — ViKj^, as well as those of aji — VjdilogN. The initial 
values of Loij = N~^doKij are given by Rij = in (33) and those of aoi = N^'^dotti by 
— {diH + HdilogN), where the simple harmonic slicing condition doN + N^H = was 
used in the latter. 

Upon inspection, one sees that all fields evolve along the do axis except those obtained 
from the wave equations and their adjoints. Dividing these equations through by N (to 
obtain proper time), just as in the case of Maxwell's theory we obtain characteristic fields 
of the form, for the choice = x^, 

■^ct[Loij ± a-\Liij - v'^Laij)], (74) 

-^a[aoi ± a~^{aii - v'^aai)]. (75) 

These fields propagate at light speed along the cone, and they have the dimensions of 
curvature, as evident from (29) and (30). At a given event, the "electric" or "tidal" parts 
of the Riemann tensor (A'"~^it!oioj) are constructed from contributions from these fields, 
which arrive from along the past cone, and from additional contributions which arrive 
from along the past timelike curve orthogonal to the slice of the given foliation in 
which the event lies. 

By writing (74) and (75) in the equivalent forms 

■^o,[N-% ± a-^(Vi - v''^a)]Kij, (76) 

■^a[N-'do ± a-\Vi - v''Va)]a^, (77) 

we see that the "characteristic fields" are differential extrinsic curvature and differential 

acceleration, with differences taken along the past cone of an event using the many pairs 
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of null generators corresponding to advanced and retarded times for all spatial directions 
in Tjt with origin at the event. 

Finally, we remark that net matter contribution to the wave equations (62) and (65) 
are included by adding to the right-hand sides the quantities 

e,, = iVr,, =iV/c[-iV-i4p,, + 2iV-iV(,p,)o + if.,(2iV'Voo + p\) (78) 

$i = Npi = -Nk [pika'' + N-^ViPoo] ■ (79) 

Note that the detailed forms of (78) and (79) depend on which Einstein equations we 
have used to display the right hand sides of (62) and (65). 



4 Higher Order Formulation 

We may ask whether a wave equation form of the Einstein evolution equations can be 
found in which both the shift and the lapse can be freely specified. The answer is yes 
pO| , but we will not give a full description here. 
Form the expression 

Aij = doQ.j + V,Vji?oo = dodoRij - 2doV{iRj)o + ViV.i^oo- (80) 

This gives a system of the form 

ndoKij = [terms with derivatives of K of order < 2 and of g of order < 3], (81) 



dogij = -2NKij, (82) 

that has no slicing term analogous to Sij. It has characteristics along the light cone and 
do- It is not quasi-diagonal because of third order derivatives of Qij arising in (81) from 
V^AA^. This system is hyperbolic non-strict |21[|. On the other hand, the fourth order 
equation for gij deduced from (81) and (82) is also not strictly hyperbolic because of the 
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appearance of Oq in its principal operator (SqD). It seems that the system, in its present 
form, cannot be written as a first order symmetrizable system; but it is potentially very 
useful, especially in perturbation theory, because the lapse and the shift need not be 
perturbed or satisfy any particular gauge conditions. The system (81) and (82) can be 
proven to be equivalent to the full Einstein equations using the twice-contracted Bianchi 
identities and the appropriate Cauchy data: Gap = and Vtij = on Sq. This system can 
be directly adapted for the study of the non-linear perturbative regime of high-frequency 
wave propagation W\ |2H . 



5 Concluding Remarks 

One of the numerous attractive reasons for writing general relativity in flux- conservative 
FOSH form, as emphasized recently in is to facilitate emulating the successful strategy 
applied in computational fluid dynamics, rather than tailoring the form of Einstein's equa- 
tions to suit each particular application. However, in fluid dynamics in flux- conservative 
FOSH form, the only characteristics that can appear are the physical ones associated 
with the "sound cone" or "flow lines," because there are no analogs of the spatial co- 
ordinate and slicing "gauges" of general relativity that give rise to constraints and can 
lead to unphysical characteristics. In this paper, by modifying the equations of motion 
of Einstein's theory in a specific way {R-ij — fijj), and by choosing time slicings in which 
the lapse function obeys a causal wave equation, we have incorporated full spatial co- 
ordinate freedom while assuring that only physical characteristics are present. Among 
other things, this means that causal boundary conditions can be imposed naturally on 
bounded regions of spacetime. Gravitational physics and its mathematical expression 
have been smoothly melded. 
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